Energy extremality in the presence of a black hole 
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l> : ABSTRACT 



We derive the so-called first law of black hole mechanics for variations about 
stationary black hole solutions to the Einstein-Maxwell equations in the absence of 
sources. That is, we prove that 5M = kSA+u8J + VdQ where the black hole parame- 



ters M, k, A, lu, J, V and Q denote mass, surface gravity, horizon area, angular velocity 
of the horizon, angular momentum, electric potential of the horizon and charge re- 



spectively. The unvaried fields are those of a stationary, charged, rotating black hole 
and the variation is to an arbitrary 'nearby' black hole which is not necessarily sta- 
tionary. Our approach is 4-dimensional in spirit and uses techniques involving Action 
variations and Noether operators. We show that the above formula holds on any 
asymptotically flat spatial 3-slice which extends from an arbitrary cross-section of 
^ ■ the (future) horizon to spatial infinity. (Thus, the existence of a bifurcation surface is 

irrelevant to our demonstration. On the other hand, the derivation assumes without 
proof that the horizon possesses at least one of the following two (related) properties: 
(i) it cannot be destroyed by arbitrarily small perturbations of the metric and other 
fields which may be present, (ii) the expansion of the null geodesic generators of the 
perturbed horizon goes to zero in the distant future.) 



I Introduction 



Working directly with the Action principle which defines a field theory typically lends 
clarity and unity to the basic formal relationships of the theory. Thus, in the case of 
symmetries and conservation laws, the first variation of the Action S, leads directly to 
the conserved "charges" themselves [p]] , while its second variation leads to conclusions 
of which the following is typical: every time-independent solution of the field equations 
is an extremum (critical-point) of the total energy E 

Two directions of generalization of this last theorem suggest themselves as nat- 
ural. On one hand, one can inquire into the behavior of the second variation, 
or "Hessian" of E (corresponding to third variation of S itself); on the other hand, 
one can allow for the presence of a spacetime boundary, such as will occur naturally 
if one narrows one's attention to the region outside a black hole horizon. Clearly, 
both directions of generalization are relevant to the stability and (more generally) 
the thermodynamics of black holes; but the latter kind of generalization seems to be 
more readily accomplished, as well as being a necessary prelude to the former kind 
in the black hole case. This second direction is pursued in the present paper, and 
leads — via steps which we now sketch in advance — directly to the so-called first law 
of black hole mechanics (or thermodynamics). 

A central concept in connection with local transformation groups is that of the 
Noether operator, an object which acts (in general as a differential operator) on 
an infinitesimal group-generator £ to give a spacetime current J a corresponding to 
the variation of S induced by £. (See equations (1) and (3) below.) When the 
group is the diffeomorphism group, £ is a vectorfield £ a , and the Noether operator 
is correspondingly a tensor (density of weight 1), with the current being formed as 
T b a ■ £ 6 . Under appropriate conditions of asymptotic flatness, a Poincare quotient- 
group of diffeomorphisms serves as a global symmetry group of the theory, and — for 
solutions of the field equations — the Noether operator yields a corresponding family 
of conserved charges when applied to the Poincare generators and integrated over an 
asymptotically flat spacelike hypersurface or "slice" .[] 

1 From a physical point of view, it seems best to interpret such a global charge, say the linear 
momentum, as the change in S due to an infinitesimal variation in which the system under consider- 
ation is translated relative to the environment whose presence is implicit in the use of asymptotically 
flat boundary conditions. Certainly, genuinely asymptotically flat metrics do not occur in nature, 
and the translation or rotation of a subsystem cannot in practice be meaningfully extended beyond 
some large but finite radius. By taking the existence of the environment more explicitly into ac- 
count, it ought to be possible to interpret the conserved charge as the change in the total S, with the 
various divergence terms one customarily adds to the Lagrangian serving only to allow a nominal 
splitting of S, and therefore of SS, into separate parts associated with the approximately isolated 
subsystem and with the environment. The same, "more realistic" approach is also important for 
correctly interpreting the divergent or conditionally convergent integrals which occur in connection 
with the angular momentum and boost conservation laws. (For the elements of such an approach 
see 0). 
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The extremality theorem cited above then follows directly from taking an ad- 
ditional variation of the defining equation (1) for the Noether operator: the inte- 
gral identity that results from the variation yields the theorem when evaluated on 
a spacetime region whose fields are chosen to interpolate between the unvaried and 
the varied configurations. In its maximal generality, this theorem asserts that any 
field- configuration which both solves the equations of motion and is invariant under 
a given infinitesimal symmetry, is an extremum of the corresponding conserved quan- 
tity; and that in fact any two of the properties, "solution" , "symmetric" , "extremum" 
imply the third.0 

What happens then if we apply the same considerations, not to the entire space- 
time manifold, but just to the region outside an event horizon? 

But first, why is it necessary to exclude the interior region at all? It is not that 
the general theorem ceases to apply just because a horizon is present, but that, for 
subtle reasons, it does not furnish the kind of information one might expect in that 
case. With the Schwarzschild metric, for example, there are two possibilities for what 
is inside the horizon. There might be the familiar second asymptotic region joined to 
the first by a "throat", or there might be no such region if the internal topology is 
such that the "throat" leads one back out to the same asymptotic region from which 
one came (the latter case resulting from the former by suitable identifications.) In the 
former case the quantity l E i which occurs in the theorem turns out to be the difference 
of the energies seen in the two asymptotic regions, rather than the effective mass of 
the black hole, while in the latter case, the meaning of E is physically appropriate, 
but the solution is no longer globally stationary technically, because the non-simple- 
connectivity of spacetime obstructs the extension of the timelike Killing vector to 
the interior region as a single-valued field. In neither case does one learn anything 
directly about variations of the physical black hole energy. However, if we truncate 
the manifold at the horizon then the remaining, external portion of the spacetime 
possesses both a single asymptotic region and a globally defined Killing vector, and 
we can expect that a suitable generalization of our theorem will furnish the kind of 
information we desire. 

Consider then, an arbitrary variation of the metric in the neighborhood of an 
asymptotically flat slice, S which meets the horizon in an arbitrary two-dimensional 
cross-section. (The freedom to choose the cross-section freely is important, because 
the ability to push such a section forward in time appears to be crucial to an under- 
standing of the second law of black hole thermodynamics; see the further remarks in 
the conclusion section below.) As long as the variation stays away from the horizon, 
the general theorem applies exactly as before, and we conclude that the energy of 

2 Similar relationships can be derived naturally via Hamilton's equations, but the Lagrangian 
approach has the advantage of leading more directly to objects with a spacetime meaning, such as 
conserved currents and fluxes (cf.Q for example). And of course, a spacetime approach is particularly 
well suited to dealing with black holes, whose very definition involves the global causal structure. 
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a time-independent black hole is an extremum against variations of the fields which 
vanish in a neighborhood of the horizon. From this alone, it follows immediately that 
for arbitrary variations (including non-stationary ones), the change in the total energy 
can depend only on the behavior of the variation on the horizon itself, or more pre- 
cisely, on its behavior near the 2-surface where the horizon meets the slice on which 
the energy is being evaluated. 

We have come this far on the basis of reasoning valid for any field theory. To 
discover specifically which horizon-variables determine 5E, we will need to use the 
concrete form of the Action S. (One might think that, in addition to horizon variables, 
the choice of spacelike hypersurface might enter into 5E, but that cannot be true, 
because the energy of the varied configuration does not depend on the slice on which 
it is evaluated — assuming that the variation itself does not lead out[] of the space 
of solutions.) For simplicity, consider the special case of a static, uncharged black 
hole in pure gravity, and let £ a be the asymptotically timelike Killing vector of the 
unvaried solution. Then, the integral for 8E that evaluated to zero in the absence of a 
boundary, becomes with a boundary present, the net flux of a fictitious energy-current 
(that of the Einstein tensor of the interpolating metric) across the horizon. This flux 
can be evaluated entirely straightforwardly using the Raychaudhuri equation, and 
yields the familiar expression kSA, where A is the horizon area and k its surface 
gravity. We thus arrive at the first law for non-rotating uncharged black holes. 

The more general case of a charged and/or rotating black hole is equally straight- 
forward, and needs no special elaboration in this introduction, except to mention 
that, in dealing with angular momentum, we will need the analog for rotation of the 
"asymptotic patching" lemma employed in 0] for the case of translation. This is 
derived in an appendix (in a version strengthened significantly over that of ||, the 
improvement being made possible by our imposition of the so-called "parity condi- 
tions" on the asymptotic fields in conjunction with the use of an improved Noether 
operator over that used in j3|). 

The only sightly tricky point in the derivation will be locating the horizon of 
the varied metric, a task which in general would demand knowledge of the entire 
future of the varied solution. In fact, what we will need for our derivation is only 
that the expansion 9 be of second order in the perturbation when evaluated on the 
correctly identified varied horizon. Although "local" in itself, this assertion refers 
implicitly to the entire future, and we do not prove it herein. Instead, we reduce it 
to either of two assumptions asserting the stability (in a certain sense) of the black 
hole horizon. Such assumptions are common in discussions of black hole dynamics, 

3 The Noether operator allows us to define a generalized energy, even for field-configurations 
which do not satisfy the Einstein equation, and we will need to use this feature in our derivation 
of SE. Thus, it may be of some interest that hypersurface independence for this more general type 
of variation still follows from the identity on which the proof is based (i.e. SE can depend on the 
choice of slice only insofar as this influences the cross-section in which the slice meets the horizon) . 
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and may perhaps be viewed as special cases of "cosmic censorship". Were they to 
fail, black hole spacetimes would have to be understood as very special types of 
mathematical solutions without relevance to reality. 

The layout of the remainder of this paper is as follows. In section 2 we recall the 
definition of the Noether operator for gravity and it's relation to conserved quantities 
I], ||, [|. The rest of the main body of the paper is devoted to the proof of the first 
law of black hole mechanics for variations from stationary 'gravitoelectric' spacetimes 
to nearby (not necessarily stationary) gravitoelectric spacetimes (by 'gravitoelectric', 
we intend gravity coupled to the source-free electromagnetic field). In section 3 we 
sketch the idea for the proof based on an 'extremum' identity and set the stage for the 
proof with some technical remarks. Section 4 gives the proof for "vacuum" solutions 
and section 5 generalizes it to the gravitoelectric case. In the concluding section, 
we discuss the significance of the result obtained and consider its relation to earlier 
work. In the appendices, the asymptotic falloffs of the metric and the electromagnetic 
potential are specified, and a lemma crucial to our proofs is proved. 

Since all variations will be about stationary solutions with vanishing spatial mo- 
mentum p, we will have 5E = 5M to first order, where M is the mass of the black 
hole (compare the relation E 2 = M 2 + p 2 ). From now on we will denote first order 
variations in energy by 5M, and often call the energy itself M. 

For a "minimalist" account of what follows, see the first part of reference [H]. 

II The Noether operator for gravity 

For the geometrical invariances of classical field theories, the 'Noetherian' relation 
between symmetries and conservation laws can be codified in terms of a differential 
operator T b which acts on an arbitrary diffeomorphism-generating vector field £ a to 
produce a conserved current T b ■ £ b . One can associate such a Noether operator to 
any first order Lagrangian density C[Q] via the identity |TJ, H 

8S = [ (5C/5Q)6Qd n x + f fT b a ■ £ b da a (1) 
Jn Jan 

In this equation, the Q are the dynamical field variables, S = j n Cd n x is the action, 
(5C/5Q) is it's variational derivative, and SS is the result of a variation in which (i) 
the variables Q 'change in place' according to 

5Q = -fCsQ (2) 

(jQ being the Lie-derivative), and (ii) the region of integration shifts its boundary 
by an amount /£ a . Notice that the arbitrary scalar field / in (1) and (2) helps to 
reduce the ambiguity of T b a , which, for / = 1, would remain a solution of (1) after 
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supplementation with an arbitrary total divergence formed from £ and the fields Q. 
It is implicit in (1) that T b = T b a [Q] depends locally on Q; and (1) is required to hold 
for arbitrary £, / and Q whether or not Q solves the equations of motion. 

The expression T b ■ £ b can be integrated to produce conserved quantities or 
'charges' associated with the asymptotic symmetries of an isolated system in a manner 
discussed in [pj, wherein the interpretation of these global charges as variations of the 
total action is also explained. Without going into further detail here, let us remark 
that the interpretation is valid when £ is an exact symmetry of the background fields, 
and the presence of the function / in (2) is what allows the variations in question to 
remain compatible with the asymptotic flatness of the metric g^ at large radii. If £ a 
generates time translation then / da a T b ■ £ fe will be the total energy; if it generates a 
rotation it will be an angular momentum component, etcetera. 

The form of S grav which leads to the Noether operator we use for gravity is the 
first order form obtained from the covariant action (1/2) / RdV via integration by 
parts. To define it in an intrinsic manner we can introduce a background connection 

o 

Va (say torsion free) and let T bc represent the difference therefrom of the metric's 

o 

connection V Q . (In a specific coordinate system, if V a is chosen as the coordinate 
derivative operator, then r£ c will just be the Christoffel symbol as usually defined.) 

The first order Action is then S grav = l/(2fc) / dV(T^ c T b ad - T a ab T b cd + R cd )g cd , where 

o o 

R c d is the Ricci tensor of Va and k = 8nG (in 4 spacetime dimensions). The Noether 
operator which answers to S grav turns out to be that acting as follows |], |J : 

T b a ■ £ fe = -G a b £ b + W h iy [a i b] + c>£ 6] ) (3) 

where T b a = (-g)~ l l 2 T b a , u a = T a bc g bc - T b bc g ac = {gy 1 Vb{-gg ab ), and x'V 1 = 
1/2 [x a y b — x b y a ). We also define 

W ab . £ C = (y[a£b] + ^6]) (3^ 

Note on Notation: In the various integrals which appear, we shall incorporate all the 
density weights into the volume elements which will accordingly be scalars (of density 
weight 0). We write the various volume elements thusly: dV denotes the element of 

For £ corresponding asymptotically to a rotation, the integral defining SS converges only con- 
ditionally and only to first order in the variation 5. This obstacle to giving angular momentum a 
precise meaning derives from the fact that S itself converges only conditionally in general. It can 
be overcome by redefining S as the Action inside a sufficiently large sphere, and taking limits as the 
size of the sphere goes to infinity 
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4-volume, dS a denotes the volume element on a 3- dimensional submanifold and dS a b 
that on a 2-dimensional submanifold (which will always be topologically a 2-sphere). 
We also absorb into the volume elements dS a , etc. the orientations which implicitly 
occur in Stokes' Theorem, although in this case the received notation does not allow 
us to indicate which orientation is being taken. In general relations like equation (9) 
we will adhere to the rule that the boundary of a region is to be oriented outward. 
In integrals over separate spacelike or null hypersurfaces like the integrals over £ in 
(13) and (14), we will orient the surfaces upward (as if they were regarded as the 
boundaries of their pasts). When the surface is a portion if of a black hole horizon, 
this has the consequence that H is oriented "inward" (as if regarded as the boundary 
of the region of spacetime outside the black hole). Other orientations will not be 
made explicit, but can be deduced from the context. 

With these conventions, the energy E of a gravitational field evaluated on a spatial 
asymptotically flat slice £ without boundary is given by 

-E = (l/k)J^T b a -edS a 

where £ is an asymptotic time translation and T b a is an ordinary tensor operator, not 
an object of density weight 1. Note that when the vacuum Einstein equations are 
satisfied, we have from (3) and (3a), 

£ = -(1/2) / W c ab -CdS ab (4) 

where <9E is the "2-sphere at spatial infinity" and k = 8tcG is set equal to 1. 
The assumed fall off conditions on the metric are given in the Appendix. 

Ill Preparation for the proof 

For a more detailed exposition of the ideas in this section see [0, [7j. 

III.l The "extremum identity" for variations about a sta- 
tionary solution, and its application to black hole space- 
times 

We notice that for a variation which corresponds to a pure diffeomorphism gener- 
ated by £ a (which corresponds to putting / = 1 in equation(l)) we have for the 
gravitational action in (3) 

5S grav := StSgrw = I dV(l/2)G ab C i9ab + f T b a - £ b dS a (5) 
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But S grav is invariant under diffeomorphisms which preserve the background connec- 



o 



tion V- Let us assume that £ preserves the background connection; then 



8 e S = 



(6) 



[ dV(l/2)G ab C ( g ab + [ 7?-£ 



b dS a = 



(7) 



Note that this identity holds for arbitrary field configurations, which do not need 
to be solutions to the field equations. We now take the variation of (7) under an 
arbitrary variation 5g ab of the metric. We get 



[ (l/2)5(dVG ab )C (9ab + I dV(l/2)G ab 5(C ( g ab ) + 5 [ T b a -i b dS a = (8) 



Let us further assume that we are varying about a metric which is both a solution 
to the field equations (G ab = 0) and stationary with respect to £ a (C^g ab = 0). Then 
the above equation reduces to 



which we shall call the extremum identity. 

For the proof, we will take g ab to be a stationary black hole solution to the vacuum 
field equations with killing vector field £ a , and we will choose f2 as the region bounded 
between the horizon '<9-B' of the stationary black hole and two asymptotically flat 
spatial slices, and Ej (Ey being to the future of Sj), each of which intersect the 
horizon at one "end" and go out to spatial infinity at the other "end". Also, we 
will choose £ a to be that Killing field of the stationary solution which is null on the 
horizon. (The general Kerr-Newman metric always admits such a Killing vector; the 
horizon "rotates rigidly".) Then 



where t a and 4> a are the Killing vectors of the stationary solution which are asymptot- 
ically generators of time translation and rotation respectively; and u is the angular 
velocity of the stationary horizon with respect to infinity. We aim to reproduce the 
first law with E, A, etc. evaluated on the slice E/. 

To that end, a suitable choice of "gauge" for Sg will be very helpful. The gauge- 
freedom arises in the following, well-known manner. If g is a metric on a manifold 




(9) 



c = 



t a + wfi 



a 



(10) 
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iV and g is a metric on N, then the pairs {N,g) and (N,~g) may be said to be 
"infinitesimally close" iff there exists a diffeomorphism $ : iV — > N, such that <&*~g = 
g + Sg where Sg is infinitesimal. Q Such a diffeomorphism $ (call it 'allowable') is not 
unique: if II is allowable and \1/ is any diffeomorphism infinitesimally close to the 
identity, IT o \|> is also allowable. We will use this freedom in the specification of $ in 
(b), (c) and (d) below. 

We now divide the 4-manifold outside OB into 3 parts: 

• Sj and its past 

• E/ and its future 

• The region Q between Ej and 

and we contemplate a preliminary choice of Sg a b as follows. 

(a) On Sj and to it's past, Sg a b = 

(b) On and to it's future (g a b + Sg a b) is a black hole solution to the linearized field 
equations with mass M+SM and angular momentum J+SJ. We assume, without loss 
of generality, that the horizons of the varied solution and of the stationary unvaried 
solution coincide in this region. 

(c) In Q: we require that the part of OB between E; and Ey be a null surface for the 
metric (g a b + Sg a b)- We shall henceforth refer to this portion of the horizon of g + 5g 
as H. 

(d) On dB it is possible to "line up" the null geodesic generators k a , for g a ^ with 
those for g a b + 5g a b, so that Sk a is proportional to k a . This alone does not imply 
that the affine parameters of the unvaried and varied null geodesic fields agree on H. 
In general they will differ by infinitesimal amounts. But it is always possible to use 
the infinitesimal diffeomorphism freedom whose existence we emphasized above, to 
identify slices in the unvaried spacetime with slices in the varied spacetime, in such a 
way that the affine parameters actually do agree. Thus, on H, we can always arrange 
for 5k a = 0, without in any way constraining the location of Ej with respect to the 
varied metric. 

Based on our choice of Sg, we shall end up proving the following statement (for 
the vacuum case): Given a stationary vacuum black hole spacetime N with metric 
g and an infinitesimally nearby pair (N,g), there exists an allowable diffeomorphism 
$o ; which maps g to g + Sg and E (a slice in the exterior of the horizon of the varied 
spacetime) to an "infinitesimally nearby" slice E such that the first law holds on E 
with regard to the infinitesimal variation Sg. But, since all slices in the stationary 
spacetime have the same horizon area, we will have actually proved the first law in the 
context of an arbitrary allowable diffeomorphism, not just $o- Similar considerations 
apply to the gravito-electric case. 

5 The notion of infinitesimal can be made precise in the usual way: replace $ and g by smooth 
parameterized families $(A) and g(X) with $(0)*p(0) = g, and define Sg := d[$(\)*g(\)]/d\\\ =0 d\. 
Similarly, we interpret (b) below (for example) to mean that the horizons coincide for the members 
of some 1-parameter family of solutions whose derivative is Sg a t 
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We choose our background connection V, so that it is compatible near spatial 
infinity with the flat metric rj a b which we will use to define our fall off conditions. 
We also choose it so that it is globally Lie derived by £. (For example we could 

o o 

use V = (1 — A)V/iat + ^stationary with A = 1 on H, A = at spatial infinity and 
£ 5 A = 0.) 

Note that the requirement that g + Sg on and to the future of X f be a solution 
with mass M + SM and angular momentum J + 5 J different in general from the mass 
M and angular momentum J of the stationary solution implies that in the region Q, 
g + Sg cannot be a solution to the field equations. 

III. 2 An important technical point: asymptotic patching 

Although we would like to make a variation satisfying the conditions (a), (b) (c) 
and (d) of section III. 1 above, we notice that such a variation would violate the 
asymptotic fall off conditions on the metric (for these conditions see the appendix). 
We can see this in the following way: The mass-information for a solution to the field 
equations is in the '1/r part' of the metric. Now, the mass of the initial and final field 
configurations on the slices Ej and £/ are M and M + 5M respectively. So in between 
these slices near spatial infinity the 1/r part has to be time dependent. But this means 
that d t (g+Sg) will fall off as 1/r and not 1/r 2 as required. Similar considerations apply 
to 'angular momentum information'. Hence, the fall off conditions in the appendix 
would have to be violated. If this happened, the variations of the action S grav would 
be ill defined, meaning that we could not use the variation outlined in section III. 1 . 
This is the technical problem. To avoid it, we employ the following trick, along the 
lines of the trick performed in 

First recall that our aim is to produce from the extremum identity the explicit 
horizon term to which we know that SM must reduce. Consider, then, an asymptot- 
ically flat spatial slice EcH extending from H to spatial infinity. For an arbitrary 
black hole solution to the field equations, the energy M evaluated on the spatial slice 
£ is given by equation (4) as M = —(1/2) J s2 W° b ■ t c dS a b (where is the 2-sphere 
at spatial infinity and t a is the asymptotic time translation vector field). 

We can extend the definition of the energy to non-solutions by setting 

- M = [ \w? ■ t c dS ab - / G a b t b dS a 

= (/, - L )\w?-t c dS ab + / l w: b -t c dS ab - [ G a b t b dS a (11) 
where S^ nH denotes the 2-sphere in which E meets H. Of course this agrees with the 
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energy on a solution because the Einstein tensor vanishes on a solution. Similarly we 
can extend the definition of the angular momentum to non-solutions by setting 

J = I V<f ■ <\> c dS ab - [ G a b ct> b dS a (12) 

(where <p a is the asymptotic rotational vector field). In view of Stokes theorem and 
eq. (3), our definitions of energy and angular momentum extended to non-solutions 
are, respectively: 

-M = [ T b a -t h dS a + [ l -Wf-t c dS ab (13) 
J* Jsl nH 2 

and 

J = [ T b a ■ (j) b dS a + I -Wf ■ (j> c dS ab (14) 

o 

Relative to a given choice of background connection V, equation (13) can be 
interpreted to say that the total energy, M, of the spacetime is the sum of the energy 
outside the horizon (the first term on the right hand side of (13)) and the energy 
of the horizon itself (the second term on the right hand side of (13)). The terms 
in equation (14) can be interpreted in a similar way. Equations (13) and (14) will 
provide the second key ingredient for our proof. 

From the lemma in the appendix and lemma 3.1 of M, we know that one can 
construct a field configuration g ab + Ag ab in a neighborhood of £/ such that: 

1. its mass as given by equation (13) is still M + 5M and it's angular momentum as 
given by equation (14) is still J + 5 J, but asymptotically it agrees exactly with g ab . 

2. near the horizon on S/ , g ab + Ag ab agrees exactly with g ab + 5g a b- 

(Notice that g + Ag then has to violate the field equations in some region far from 
the horizon.) 

Using these facts we can construct the new variation Ag in Q under which we 
choose to evaluate equation (9). Ag is such that: 

[i] Ag = at Si 

[ii] Ag agrees with 1. and 2. above; i.e. g + Ag has mass M + 5M and angular 
momentum J + 5 J on Ej, and Ag near vanishes near spatial infinity. 

[iii] Ag within Q vanishes in a neighborhood of spatial infinity. In other words Ag is 
of compact support within Q. Also, H remains a null surface for the metric g ab + Ag ab , 
and Ag ab = 5g ab near H. 

From [iii] just above and (d) of section III. 1 , we have that on H 

Ak a = (15) 
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From this property of the variation in conjunction with its other properties one ob- 
serves the following two useful facts: Since H is still a null surface its null normals 
must be proportional to g a bk b — k a . So 

Ak a = ck a (16) 
for some real c. And since £ is null on H, 

C = ak a (17) 
for some coefficient of proportionality a. 

To summarize: The trick is to replace 5g by Ag. This allows the variation to be of 
compact support, freeing us from having to evaluate variations near spatial infinity 
As a final remark, we note that the variations 5 and A are identical near H and 
we shall continue to denote by 5, the variations of quantities near H. 

IV The proof for gravity alone 

We would like to emphasize that what we have described in the preceding sections 
is just the necessary background to do the actual proof of the first law. Once this 
background is assumed and absorbed, the proof itself is extremely simple (cf. 0). 

The extremum identity, equation (9), written for the A variation is 
A Jan ■ £, b dS a = 0. But Ag = near spatial infinity and on Sj, whence 

A / T b a -edS a + 5 f T b a -edS a = (18) 
Taking the variation of equations (13) and (14) with S = Sj we get 

— 5M = A [ T b a -t b dS a + 5 I l W : h -t c dS ab (19) 



and 



5 J = A [ T b a - <p b dS a + 5 [ \w? ■ <p c dS ab (20) 

J *f J %nH 2 



Equations (10), (18), (19) and (20) yield then 

- SM + u8J = 8 [ \wf ■ edS ab - 5 [ T b a - i b dS a (21) 
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Using the definition of the Noether operator for gravity and Stokes theorem (and 
recalling that Sg = on Ej H H) we get 



SM + ujSJ = 5 Gi£ b dS a (22) 

H 



a relation of some interest in its own right. Notice that (22) can be read as equating the 
change in energy or angular momentum to what would be the corresponding "energy 
flux" or "angular momentum flux" through the horizon if SG a b were interpreted as 
the stress energy of some fictitious matter field. 

We also note that this relation shows in another way that the left hand side of 
(22) must reduce to a surface integral on E/ D H, because on the right hand side we 
can confine S to a neighborhood of E/ D H on H, as in [1], and in the limit can put 
5g a b oc 9(X — A/), where 6 is the step function. Yet another way to show this, perhaps 
the most systematic of all, would be to use the theorem in J5J to derive a "potential" 
for 5{Gl^ b ) from the identity d a 8( y f^gGl£ s b ) = 0, see the discussion in ||. 

In evaluating the right hand side of (22), let A be an affine parameter for the null 
geodesic generators of H so that k a = dx a /d\. Then 

dS a = -k a d\d 2 A (23) 

where d 2 A is the area element of the 2-sphere cross sections of H (and A increases 
toward the future). Moreover, by differentiating equation (17), it follows immediately 
(in view of the fact that k a is the tangent to an affinely parameterized null geodesic, 
whence k a V a k b = 0) that 

da/dX = k (24) 
where k is the surface gravity of the stationary black hole, defined by 

e^ a e = (25) 

In the extremum identity, £ is not to be varied, i.e. A£ = 0. Then since 5k a = 
as well, the right hand side of (22) becomes 



ak a k b {8{G ah d 2 A))dX (26) 

H 



Using G a b = for the stationary solution and the fact that k a remains null under the 
variation, we get 

- fd 2 Aj dX ak a k b {5R ab ) (27) 
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To evaluate (5R ab )k a k b we use the Raychaudhuri equation for the rate of expansion 
of a congruence of null geodesies (see for example |§): 

dO/dX = -(l/2)6 2 - a ab a ab + u ab u ab - R ab k a k b (28) 

where 9 is the expansion, a the shear and u> the twist of the null geodesic congruence 
with k a parameterized by affine parameter A. We apply this equation to the null 
geodesic congruence on H (for which anyway u = 0) and take its 5 variation. Using 
the fact that the expansion and shear of the stationary solution vanish we get: 

d{56)/d\ = -{5R ab )k a k b (29) 

Inserting this into (27) yields 



5 J G^ b dS a = j d 2 A J d\ad(56)/d\ (30) 

= (fd 2 Aa56)\{-fd 2 Ajd\^-56 (31) 

where i,f refer to the 2-spheres which are the intersections of H with E$ and 
respectively. 

But 59i = by construction, and 88 f = under the assumptions of Lemma 1 
(which is given at the end of this section). Also, 9 = (l/d 2 A) d(d 2 A)/d\ by definition. 
Hence, 



5 



Using (|24|) together with the fact that for the stationary solution d(d 2 A)/d\ = 0, we 
get for (32) 



* A I ""-h'^sr (33) 



But 5X = 0, so 



5 J H G a b i b dS a = ^fJdX 



d(5d 2 A) 



dX 



K8{jd 2 A)\l 



K5Ud 2 A)\ f 



(34) 
(35) 
(36) 
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where in the last line we have used that the variation 5 vanishes on Ej. Hence we get 
by (22) 

8M = k5A + u5J (37) 

where everything is evaluated on Ej. 

To complete our proof, we have to justify setting 89 f = 0, which we do via the 
following lemma. 

Lemma 1: Let (g a b, {0}) be a stationary black hole solution to the field equations, 
where {</>} denotes all fields other than the gravitational field, and let the stress-energy 
tensor T ab of {0} satisfy the positivity condition that T ab k a kb > for every null k a . 

Consider the spacetime region Q' to the future of some spatial slice E/. Denote the 
part of the horizon in this region by H'. Let [g a b + 5g a b, {<p + 5<j)}) be an infinitesimally 
differing solution, with the relevant part of its horizon coinciding with H' . 

Then the variation 58 of the expansion evaluated on Ej vanishes if either of the 
following (rather weak) stability conditions holds: 

(i) It is not the case that arbitrarily small perturbations of the metric and the 
other fields can destroy the horizon; or 

(ii) 59 approaches zero in the distant future along the horizon. 
Proof: 

(i) (9ab + bgab) {0 + <50}) is a solution to the linearized field equations. Therefore 
(g a b — 8g a b, {0 — 5(p}) is also a solution to the linearized equations. Thus, if for the 
former solution, the variation of the expansion on H' is 59 then for the latter solution, 
the variation in the expansion will be —59. 

Since 9 = for the stationary solution, there exists a solution to the full Einstein 
equations in a neighborhood of the stationary solution whose expansion at first order 
is 59, which, as just explained, can be arranged to be negative unless it is always 
zero. But it is known that a negative expansion implies a singular evolution of the 
horizon. Namely, by using ("to all orders") the Raychaudhuri equation (28), the 
Einstein equation and the positive energy condition above, one can conclude that if 
the expansion 9 = 9 is negative at some A = A (the geodesies are converging), then 
a conjugate point (caustic) will develop at or before some definite future instant Ai, 
where Ai depends on 9q. But the occurrence of a conjugate point contradicts the 
definition of the horizon, because no geodesic generator can remain on the horizon 
after reaching a conjugate point. Thus, the only possibility is that the generator in 
question terminates in a singularity, i.e. that the perturbation destroys the horizon, 
contradicting assumption (i). 

(ii) Applying the variation 5 to the Raychaudhuri equation and noticing as before 
that the expansion, shear and twist of the stationary solution vanish, we get from 

<m 

d(59)/d\ = -{5R ab )k a k b = -5T ab k a k b , (38) 
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in view of the Einstein equation. But from (28) and the Einstein equations it follows 
that for the stationary solution 

R ab k a k b = T ab k a k b = 
Hence the positivity condition on the stress energy tensor implies that 

5T ab k a k b > 

for arbitrary variations of the fields. But as 6T ab is linear in the infinitesimal variations 
(6g ab , {50}), this is impossible unless 6T ab k a k b = for all variations, whence 

d(86)/d\ = (39) 

for all variations. Hence 89 f = 59^ where £#00 is 59 evaluated at large affine parameter 
and 69 f is 69 evaluated at £/. But 59^ = by assumption. Hence 69 f vanishes as 
well. (Notice that if £ f intersected the horizon in a bifurcation surface as in || then, 
since by definition a = at such a surface, a59f = would be immediate and lemma 
1 would not be needed.) 

V The proof for gravito-electric spacetimes 
V.l The Electromagnetic Noether operator 

In order to generalize (37) to gravito-electric spacetimes, we will apply the extremum 
identity (9) to field variations about a stationary, charged black hole solution of the 
Einstein-Maxwell equations. From the general argument in section 1, we know that 
8M and 6 J will be expressible entirely in terms of "horizon variables" , and since the 
total Action is now the sum of Sg Tav cincl a single electromagnetic term, we know that 
a single extra term will appear in (37). We can find this term by using the gravito- 
electric versions of equations (9), (13) and (14). To do this, we need expressions for 
the relevant Noether operators. We display them below (for details see Q). 

The Action for the electromagnetic field is S em = (—1/4) / F ab F ab dV where F ab = 
V ' a A b — V b A a is the electromagnetic field tensor and A b its potential. The Action for 
the combined gravitational and electromagnetic fields is S grav + S em (S gra v has been 
defined in section II). The action of the gravito-electric Noether operator on a vector 
field £ fe is: 

1~total b ' ^ := V~ 9 ^ total b ' 6 = V~ 9 (^Tm b + ^grav b) ' £ & (40) 

where Tg rav b ■ £ b is given in equation (3) and T em l is given |§ by 

TL b -e = F ab ^A b -{l/A)F cd F cd C (41) 
= T a Sb e-V b {i c A c )F ab , (42) 
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Tg b being the stress energy tensor of the electromagnetic field, 

T ab = F ac F b ^ _ ^ 1 g ^F cd F cd (43) 

The electromagnetic analogue of (3) is therefore 

T a ernb ■ e = {T a Sb + {V c F ac )A b )e - V b {F ah A c C) (44) 

and the total Noether current is — as in (3) — the sum of a divergence with a term 
which vanishes "on shell": 

T t a otalb -e = [T a Sb -G a b + {V c F^)A b \e 

+V b (Wf ■ C - F ab A c C) (45) 

wiihWf-£ c given by (3). 

V.2 The extremum identity 

The extremum identity reads now 

* fT t a otalb -edS a = (46) 

JdQ 

As before, we will apply this with £ a = t a + u<f) a where t a and <p a are Killing fields 
which are asymptotically time translational and rotational respectively; £ a is that 
killing field which is null on the horizon and uj is the angular velocity of the black 
hole. 

V.3 Some useful properties of the stationary solution 

The variations will be made around a stationary black hole configuration, g ab , A c . 
The following properties of this stationary black hole solution will be of use: 

(a) C^A a = 0. (This amounts to a suitable gauge choice and was implicitly 
assumed when we wrote down the extremum identity above) 

(b) Claim: On the horizon F ab £ a = 7^ fc for some function 7 . 

Proof: Using (28) for the stationary solution on the horizon, the Einstein equations 
and the proportionality of the null geodesic horizon-generators to the Killing field £, 
one gets 

T Sab ^ b = 

on the horizon. Substituting the explicit form of Ts ab from (43) into this equation, one 
concludes that F ac C, c is a null vector. From the antisymmetry of F ac , one concludes 
that this null vector has to be proportional to £ a , since it is both null and orthogonal 
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to f. 

(c) The 'electric potential' —^ a A a is constant on the horizon. (Contract the 
identity ^ a F a b = C^Ab — Vb((, a A a ) with any v b tangent to the horizon and use (a).) 
Following the common notation, we will write V = —^ a A a for this constant. 

These relationships express the fact that the horizon behaves like a conducting 
surface : in equilibrium (condition a) the electric potential must be constant (condi- 
tion c) and hence the "electric field" must be perpendicular to the surface (condition 
b) (see, for example jlQfl). 

V.4 The variation 

We use the same notation as in previous sections for spacetime regions, volume ele- 
ments and anything involving only the metric. 

The conditions on variations' are similar to those in section III. 1 . We choose 
(SA C , Sg a b) as follows: 

(a) On Ej and to it's past, 5g a b = 5A C = (we don't want to involve £j at all) 

(b) On £/ and to it's future (g a b + Sg a b, A c + 5A C ) is a black hole solution to the field 
equations with energy M + SM angular momentum J + 5 J and charge Q + 5Q. We 
arrange that the horizons of the varied and unvaried solutions coincide. 

(c) , (d) These conditions are identical to conditions (c) and (d) of Section III. 1. 

o 

We also choose our background connection, V, as in III. 1 . As before, we have to 
deal with non-solutions in the region between £, and Sj. 

V.5 The variation 'A' 

We face a similar technical problem as in Section III. 2. To get around it, we repeat 
the steps followed there. We first define expressions for the energy M and angular 
momentum J of a field configuration on an asymptotically flat spatial slice £ which 
reduce to the correct values when the source free Einstein-Maxwell equations are 
satisfied, but which are also useful for fields which are not solutions (these expressions 
follow directly from (|45| ) and the fact that the total energy and angular momentum 
on a boundary-less slice can always be written as J s T b otal a ■ £, a dSb with £ the relevant 
generator), to wit: 

-M= f T t b otal a ■ t a dS b + I l -F ab AXdS ab + / \w? ■ t c dS ab (47) 
J = f T b total a ■ rdS b + I \F ab A c <p c dS ab +<f \wf- <p c dS ab (48) 
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Using arguments similar to those in section III. 2 one can introduce variations of 
compact support of all the fields (i.e. A and g) such that the new field configuration 
on Ej has mass M + 5M and angular momentum J + 5 J, and is identical with 
(A + 5A, g + 5g) near the horizon. We will not display the relevant patching arguments 
since they are similar to those of Appendix B. We simply note that in order to 
complete the patching arguments, one has to make use of the asymptotic conditions 
on A c as well. These conditions are listed in Appendix C. 

Let us, as before, call this new variation of compact support the A variation. 
Just as earlier, it satisfies the analogs of properties [i]-[iii] of Section III. 2, as well as 
equations (15)-(17). And, we can again continue to denote A on the horizon by 5, as 
there is no difference between the two variations there. 

V.6 The proof 

We have the extremum identity: 

A f T t a otal b ■ ?dS a = (49) 
Jan 

=► A / T t a otal b ■ ?dS a + 5 I T t a otal b ■ edS a = (50) 

Using the above expressions for M and J with S = E/, and taking respectively £ = t 
and £ = in fl50D yields, in view of (f4~0"|), 



-5M = 5 ( \w? ■ t c dS ab - 6 f T a grm b ■ t b dS a - 6 f T a emb - t b dS a 



+ 6 / -F ab A c t c dS a b (51) 



5 J = 5 f \w? ■ <j> c dS ab - 5 [ T« rav b ■ 4> b dS a - 5 [ T: mb - <j> b dS a 



+ 6 / -F^A^dSat (52) 

JS r .r-.tr ^ 



Then the same manipulations of the gravitational terms in these equations as we 
performed in section IV produce here 

- 5M + u5J = (grav) + (em) (53) 
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where 



(grav) = 5 J G a b £ b dS a = -k5A + £ d 2 Aa56 



(54) 



is the same expression as before (with £ a again denoting t a + w<f) a ), and 

(em) = 5(- [ T e a m b ■ ?dS a + \i F ab A c CdS ab ) (55) 
Jh 2 Js 

Here A denotes the area of the 2-sphere S : = T,ff]H and j> s denotes integration over 
this 2-sphere. Notice that in (|55"D it is not convenient to convert the integral over 
X f fl H to an integral over H as we did with the corresponding integral for gravity in 
equation (21). 

As expected, our expression for — d~M + uSJ is just the same as before, with the 
exception of the additional term '(em)' in equation (0). To convert this term to a 
surface integral over S (as we already know must be possible on general grounds), we 
proceed as follows. Direct substitution of eq. ( f4"2"D into (55) yields 

(em) = h + I 2 + h (56) 

where 

h :=-5f H F ab CsA b dS a 
h :=5J H \F cd F cd ^dS a 
h :=5§ s A c i c \F ab dS ab . 

Now ^3 is already a surface integral over S, and I2 is easily seen to vanish in light of 
(23), (16) and the fact that £ • £ and 5£ a both vanish. For Ii, we can transform its 
integrand using (V3.a) to get 

F ab C^5A b = C^F ab 5A b ) = V ' c Y ac - V c {F ab 5A b )C 

where Y ac = F ab 5A b ^ c — (a <-> c). The final term does not contribute to the integral 
Ji because ^ a dS a = by (17) and (23), and the first term yields an integral which 
can be converted by Stokes' theorem to obtain 

h = - f s F ab 5A b CdS ac 

(remember that 5 = on Ej). Returning to ^3 for a moment, we can write 

I 3 = A C CS J l -F ab dS ab + J l -dS ab F ab 5(A c C) (57) 
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where we have used (V.3c) to take —V — A - £ out from under the integral sign. It 
is then a matter of a few lines of straightforward algebra to show that the second 
integral in ( |5TD cancels with Ji. (The algebra uses (23), (V.3b), and the fact that the 
surface element dS ab of S can be written as 

dS ab = d 2 A (kj b - k b l a ) (58) 

where l a is any future-null vector orthogonal to 5* such that k a l a = — 1.) Thus we are 
left with 

(em) = h + h + h = A£ c 8 I -F ab dS ab 

Js 2 

or 

(em) = - V5Q (59) 
because of the perfectly general relation 

Q=i l -F ab dS ab (60) 
Js 2 

Q being the charge of the black hole, evaluated as an (outward) flux through S = 
H H £/. 

Substituting this back into equation fl53|) , we get 

-5M + uj5J = -k5A + £ d 2 Aa56 - V5Q (61) 
Finally, Lemma 1 comes to our aid to tell us once again that 58 = on S, so, 



putting this back in equation (61), we get our generalized extremality theorem for 
the gravito-electric case: 

5M = k5A + uj5J + V5Q (62) 



VI Conclusions 

As we intimated in the Introduction, the primary impetus for writing this paper 
came from our curiosity about how the extremality theorem of reference |Q would 
be modified by adaptation to the spacetime region exterior to a black hole. Here 
we would like to comment more at length on the wider significance of the results 
obtained, and more generally on that of the first law itself. In the course of these 
comments we also will mention some further related work which either has been done 
already, or would seem to be worth doing. 
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As its name makes clear, the most important implications of the so-called first law 
pertain to the thermodynamic attributes of black holes; however the associated issue 
of stability is important in its own right (and arises already in the purely classical 
setting). 

By virtue of being time-independent, a stationary black hole solution may be char- 
acterized as a state of equilibrium of a self-gravitating system. If this equilibrium is 
to be stable in the thermodynamic sense, then of course it must (at least locally) be 
a maximum of the entropy S (we use 'S" for entropy in this section). Conversely, if 
some configuration does maximize S, then the second law of thermodynamics implies 
that the equilibrium state in question is indeed a stable one. Assuming that S is 
a smooth function on the relevant space of "configurations" or "states", a more or 
less necessary and sufficient condition for equilibrium is thus that S have vanishing 
gradient and negative definite Hessian. More generally, one may just define an equi- 
librium configuration to be an extremum of S, and diagnose its (thermodynamic or 
"secular") stability by the behavior of the Hessian there. 

Actually, the criteria we have just stated need to be qualified, as equilibria need not 
be unconditional maxima of S, but only so at fixed values of the relevant conserved 
quantities. In the black hole situation, the relevant quantities will usually be the 
energy E, the angular momentum J, and the electric charge Q. If we identify the 
area A with the entropy (up to a numerical factor), then a necessary condition for 
a black hole configuration to represent a thermodynamic equilibrium state is that a 
relation of the form (62) (i.e. dA = f3dE — f3udJ — f3VdQ) hold for arbitrary variations 
of the fields. (Clearly this equation entails extremality of S at fixed E, J and Q; and 
conversely, extremality implies a relation of this general form, although of course, it 
does not give us the specific expressions (24), (V.3c), etc. for the coefficients /3, V, u, 
etc. which appear in (62).) This relation to thermodynamic equilibrium is the main 
reason why it is important that (62) hold more generally than just for variations from 
one stationary solution to another. 

Although the satisfaction of equation (62) qualifies a black hole solution as an 
equilibrium state in the above sense, it has nothing directly to say about stability. 
As we just pointed out, this concerns the second derivatives of the entropy, or in 
the classical limit, of the horizon area. In thus using A as the measure of classical 
stability, we of course presuppose the more general identification of area with entropy, 
which rests first of all on (62) itself. In the h — > limit, this geometrical contribu- 
tion dominates all other sources of entropy, whence the latter can be ignored. The 
consequent fact that A can be used to diagnose classical stability is consistent with 
the interpretation of the classical law of area increase as expressing the second law of 
thermodynamics for black holes, but unfortunately the area law can be justified as an 
independent classical fact only to the extent that "cosmic censorship" holds — one of 
many examples of the close interconnection between cosmic censorship and black hole 
thermodynamics. If we could establish cosmic censorship, or prove the area law on 
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some independent basis then area maximization could be used as a criterion of stabil- 
ity independently of any thermodynamic argument: all that matters is that horizon 
area be (classically) non- decreasing, and therefore a valid "Lyapunov functional" . 

In any case, in order to investigate stability on the basis of horizon area, one should 
determine the negative-defmiteness (or lack thereof) of the "reduced Hessian" , d 2 A — 
(3d 2 E, or more generally of d 2 A — f3d 2 E — (3Vd 2 Q — (3ud 2 J. One might hope to prove 
negative-defmiteness in all directions except those corresponding to "super-radiant 
modes" in the rotating case. We have not studied this Hessian directly; however there 
exists an indirect way to investigate its eigensigns, namely the so-called turning point 
method commonly employed in studies of stellar stability and the stability of other 
astrophysical objects (cf.JTl], [T^, |13|]). This method can prove instability, but can only 
make stability plausible. As applied to the Reissner-Nordstrom and Kerr-Newman 
families of black holes, it indicates stability, or rather it indicates that nowhere along 
the sequence does a new unstable mode come into being and remain for a finite range 
of parameter value ||12|| . Indirectly, this implies that the above reduced Hessian is 
indeed negative in the required sense. A direct confirmation would be of considerable 
interest. 

We believe that our method of proving the "extended first law" for Einstein gravity 
is simpler than previous ones (cf. 0), and that for this reason, it clarifies as well 
the origin of the "unextended first law" , the one applying only to variations from one 
stationary solution to another. Aside from this possibly greater simplicity and the 
associated thoroughgoing spacetime character of the definitions and the derivation, 
the most important way in which our result differs technically from earlier ones is 
that (without making any reference to a possible "bifurcation surface" ) it allows the 
hypersurface £ on which the energy, area, etc. are evaluated to intersect the horizon 
in an arbitrary manner, and in particular to intersect it arbitrarily far into the future. 
(The importance of this type of hypersurface freedom was stressed also in |14|]-) This 
freedom in the choice of E is mandated physically, since it is only some future-segment 
of a Schwarzschild metric (say) that is relevant to an astrophysically realistic black 
hole like one formed in stellar collapse. But the more important issue of principle, in 
our view, has to do with the second law itself. 

For the theory of black hole thermodynamics to achieve a satisfactory status, it 
will be necessary not only to derive the equilibrium value of the entropy from first 
principles, but also to explain why the law of entropy increase continues to hold when 
black holes are present. One attempt to foresee how such a proof would go [[0|, 
presupposes that (say in a setting which is near enough to classical that an approxi- 
mately well defined identification of individual spacelike hypersurfaces can be made) 
the total entropy in question is that of an effective quantum density-operator asso- 
ciated to the portion of each such hypersurface lying outside the black hole(s). The 
second law is then identified with the assertion that this entropy is weakly increas- 
ing as the hypersurface advances in time along the horizon (or at spatial infinity). 
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For consistency, then, it would be necessary that our earlier considerations on area 
extremization hold for such hypersurfaces, as well as for ones tied, for example, to a 
bifurcation 2-surface. In fact, the desire to prove extremality for such hypersurfaces 
was our second main reason for writing this paper. 

Extremal theorems can also be approached from a Hamiltonian point of view, as 
we mentioned in the Introduction. Recently, such an approach was employed in ||, 
the result being a theorem similar to ours, except that the hypersurface on which the 
canonical variables were defined was assumed to meet the horizon in a bifurcation 



2-surface. Other recent papers dealing with these same general issues are [[L6| [L7|, [Lj| , 
the last two of which are rather similar in spirit to our own. Some of these are 
concerned with higher derivative or higher dimensional gravity. There is little about 
our methods which would not seem to apply in those cases as well, but that is a 
question we have not looked into. 

In concluding, let us mention a final technical point and a final issue for further 
study. In our derivation we have defined the energy and angular momentum directly 
from the Noether operator expressions described in Section II. For completeness, one 
would like to confirm that these definitions agree in general with the ADM ones, given 
appropriate asymptotic falloffs. This has been done and will be reported by one of 
us in another place ]19[. Related to this way of introducing conserved quantities, is 



the fact that E, for example, splits naturally into a sum of an "exterior energy" — a 
spatial volume integral — and a surface contribution from the horizon (see equations 
(13) and (47), or (14) and (48) for J). This raises the question whether (perhaps via a 

o 

suitable choice of background connection V) these separate terms (or their variations) 
might have any individual significance, so that one could in this way meaningfully 
disentangle the energy or angular momentum "of the black hole" from that "of the 
exterior matter". 

The work of RDS was partially supported by NSF grant PHY 9307570. 

Appendix 



A The asymptotic fall off conditions on the metric 

The following notations/conventions will be used in the rest of the Appendix. Fix 
an asymptotically Minkowskian coordinate system (t, x, y, z) in a neighborhood of 
spatial infinity. Greek indices will denote 4-d spacetime Cartesian components, Latin 
indices i, j, k... will denote 3-d spatial Cartesian components and Latin letters a, b, c... 
will be abstract indices. The matrix rj^ representing the fixed flat metric at spatial 
infinity will be diagonal with components (—1,1,1,1). We also choose our fixed back- 

o 

ground connection V M to coincide with the Cartesian coordinate derivative <9 M in a 
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neighborhood of spatial infinity. 

We define the inversion map, X, on a neighborhood of spatial infinity in the fol- 
lowing way: 

X(aJ») = -a? (63) 

The point-map X naturally induces a map on tensor fields in a neighborhood of spatial 
infinity which we shall also call X. Our falloff conditions will entail that, to leading 
nontrivial order, 

Xg a f> = g a b 



XVcSU = VcQab (64) 

The full set fall off conditions on the spacetime metric g a b at spatial infinity in 
terms of r = (x 2 + y 2 + z 2 ) 1 ^ 2 are as follows (spatial infinity being approached as 
r — ► oo at fixed t). 
Let := - 77^. Then 

(a) V = a^{x T )/r + 0(l/r 2 ), 

where a llv {x T ) is bounded and a^ v {— x T ) = a IMU (x T ) 

(b) a V = /W* T )/r 2 + O(l/r 3 ), 

where /3 aflv (x T ) is bounded and j3 a ^ v {—x T ) = —/3 aflu (x T ) 

(c) «9 Q ^V = 0(l/r 3 ) 

Now consider an asymptotically flat spatial slice, E, which is asymptotically a 
t = a slice (a being some constant). Conditions (a) and (b) above relate fields on 
the t = a slice to those on the t = —a slice, but for our estimates, we would like 
conditions dealing exclusively with fields on a single slice. The fields at (x\t = —a) 
are related to those at (x l ,t = a) by an integration of their derivatives over a finite 
time interval (At = 2a). Using this fact we obtain, on a fixed slice S: 
(al) hp, = cvCx'Vr + 0(l/r 2 ) 
(bl) <9 a V = Pa^)/r 2 + 0(l/r 3 ) 

where a^ u and /3 afll/ in (a) and (b) have been restricted to E and are respectively even 
and odd functions of x l for large r. 

The evenness conditions expressed in fl64D , (a), (b), (al) and (bl) are often called 
"parity conditions". In this appendix, we will employ the term more generally to 
denote our falloff conditions as a whole. 



B Asymptotic patching and angular momentum 

Lemma: Let g a b and g a b be two metrics which satisfy the parity conditions of part A 
of this appendix, and set J(g a b) '■= / s T a 6 - 4> a dSb where T b a is the gravitational Noether 
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operator and (j) a is an asymptotically Killing, rotational vector field which commutes 



o 



with V in a neighborhood of infinity. In a neighborhood of E, let 



9ab ■ 



bg a b + (1 - b)g ab 



where the "patching function" b(x) is defined as follows. Fix a smooth function / 
on [0, oo] such that f(y) = 1 for y < 1 and f(y) = for y > 2. Choose a radius R 
large enough so that the asymptotic coordinate system is defined for r > R/2, and 
set b(x) = f(r/R). (Strictly speaking, r might not be defined throughout the region 
enclosed by the r = R/2 two-sphere, but we still put b = 1 on that region.) Then the 
difference J(g ab ) — J(g a b) can be made as small as desired by choosing R big enough. 

Proof: A useful expression for T£ ■ M is: 



T; • r = --(V v cj> p )[2Tl a g^g°» - Vj^Tf + g^g Tp ) + K a g^g^ - I^Vl 



(In this equation, 'i?' stands for the Ricci scalar, of course; it is not the radius 
parameter involved in the definition of the patching function.) Since the second 
derivatives ddg cancel out in the second term of (]65D , the main term to worry about 
is the first, which might in principle decay as slowly as 0(l/r 2 ) because dip is 0(1) 
and T is 0(l/r 2 ). Thus the crucial fact about the first term is that the coefficient of 
Vjj^p is symmetric in 77 and p. Using this, and taking R big enough so that 0°" is a 
Killing vector for the background flat metric rj^ for r > R one has symbolically (in 
the spirit of 0]), 



As a result, the contributions to / are of the following type (as could also have been 



(Note: In the above all indices have been suppressed and by h we mean either h or 



1 



(65) 



I = J(g)-J(g)= dgdgcj) + / d<j>(g-Tj)dg 

Jr>R Jr>R 



(66) 




h) 
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Now since db = 0(1/ r), none of these integrals can be worse than logarithmically 
divergent. In fact we now show that they all are bounded and vanish as R — > oo. We 
use heavily that b and the ot^ v are even functions on the 2-sphere at constant r and 
that the d^b, j3 afll/ and M are odd functions on the 2-sphere at constant r. We denote 
the element of solid angle on the 2-sphere by dQ. Consider l\ as a typical example. 
We have 



h 



2R 



bbdhdhi 



R 



f 2R r 2R 

< / bb((3/r 2 )(f3/r 2 )<f)r 2 drdtt + / CbbCJ^U^drdQ 

JR JR 



where for ease of writing we have omitted absolute value signs around the the members 
of the inequality. Using the parity conditions for b,/3,(p and that |0^| < 2R, \b\ < 1 
we get: 

r-2R 

h < / (odd f unction /r 3 )r 2 drdtt + 8ttCJR (67) 

JR 

I\ < A\/R for some constant A\. In a similar fashion, for i = 2 to 5 are all 
bounded by some A^/R. (To show this, one also has to use, in addition to parity 
arguments, the fact that \db\ < constant /R, which follows directly from its definition 
as f(r/R).) Using parity arguments it is similarly easy to show that Iq, are also 
bounded by some Aq/R and A7/R respectively. Therefore, as R — » 00 all the above 
contributions to / vanish, and we have proved the lemma. 

Remark: In relation to the case of a slice without boundary, we could also apply 
this lemma to a proof of an extremum theorem for angular momentum along the lines 
of the extremum theorem for mass in (cf. ||). Notice in this connection that, 
our asymptotic conditions are about as weak as one could expect to be physically 
relevant, and as such are less restrictive than those used in the proof of the extremum 
theorem for angular momentum given in |4j]. 

C The asymptotic fall off conditions on A a 

We use the notation of section A above. As there, our fall off conditions will involve 
parity conditions under the inversion map X, requiring to leading order in (1/r) that 

lA a = -A a (68) 

More specifically our conditions are 
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(a) = Xn( xT )/ r + 0{l/r 2 ) where Xn( xT ) is bounded and x T ) = Xn( xT ) 

(b) d v Ap = p^{x T )/r 2 + 0(l/r 3 ) where p^ is bounded and p llu {-x T ) = -p llu {x T ) 

(c) d a d p A, = 0(l/r 3 ) 

The same reasoning as in A shows that on a l t =constant' slice E, (a) and (b) 
imply: 

(al) = xM)l r + 0{l/r 2 ) 
(bl) d^ = p^)/r 2 + 0(l/r 3 ) 

where Xn an d P^v have been restricted to E and are respectively even and odd func- 
tions on the r =constant, t =constant 2-sphere for large r. 

Finally, we may observe that all the arguments in the paper would go through if 
we imposed the opposite parity on A a from that expressed in (68). Our choice of sign 
was made only to conform with the character of the field due to a static point charge 
at the origin of Minkowski spacetime. It means in effect that the spacetime inversion 
J has been chosen to act as PT rather than as CPT. 
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